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The dynamics of a propagating magnetic island interacting with a resis- 
tive wah or a static external resonant magnetic perturbation is investigated 



C/3 

r~| ■ using two-fluid, drift-MHD (magnetohydrodynamical) theory in slab geome- 

try. In both cases, the island equation of motion is found to take exactly the 



^ . same form as that predicted by single-fluid MHD theory. Three separate ion 

polarization terms are found in the Rutherford island width evolution equa- 



(•~^ . tion. The flrst is the drift-MHD polarization term for an isolated island, and 

^D . is completely unaffected by the interaction with a wall or magnetic perturba- 

^^ . tion. Next, there is the polarization term due to interaction with a wall or 



magnetic perturbation which is predicted by single-fluid MHD theory. This 
term is always destabilizing. Finally, there is a hybrid of the other two po- 



^ . larization terms. The sign of this term depends on many factors. However, 

under normal circumstances, it is stabilizing if the unperturbed island prop- 
agates in the ion diamagnetic direction (in the lab. frame), and destabilizing 
if it propagates in the electron diamagnetic direction. 
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I. INTRODUCTION 

Tearing modes are magnetohydrodynamical (MHD) instabilities which often hmit fusion 
plasma performance in magnetic confinement devices relying on nested toroidal magnetic 
flux-surfaces.^ As the name suggests, "tearing" modes tear and reconnect magnetic field- 
lines, in the process converting nested toroidal flux-surfaces into helical magnetic islands. 
Such islands degrade plasma conflnement because heat and particles are able to travel ra- 
dially from one side of an island to another by flowing along magnetic fleld-lines, which is a 
relatively fast process, instead of having to diffuse across magnetic flux-surfaces, which is a 
relatively slow process.^ 

The interaction of rotating magnetic islands with resistive walls^^^^ and external res- 
onant magnetic perturbations^'^'^^^^^ has been the subject of a great deal of research in 
the magnetic fusion community, since such interactions can have a highly deleterious ef- 
fect on plasma conflnement. This paper focuses on the ion polarization corrections to the 
Rutherford island width evolution equation^^ which arise from the highly sheared ion flow 
profiles generated around magnetic islands whose rotation frequencies are shifted by interac- 
tion with either resistive walls or external magnetic perturbations. According to single-fluid 
MHD (magnetohydrodynamical) theory,^'^^ such polarization corrections are always desta- 
bilizing. The aim of this paper is to evaluate the ion polarization corrections using two-fluid, 
drift-MHD theory, which is far more relevant to present-day magnetic confinement devices 
than single-fiuid theory. This goal is achieved by extending the analysis of the companion 
paper, ^^ which investigates the dynamics of an isolated m.a.gneiic island in slab geometry us- 
ing two-fiuid, drift-MHD theory. For the sake of simplicity, we shall restrict our investigation 
to slab geometry. 



II. REDUCED EQUATIONS 
A. Basic equations 

Standard right-handed Cartesian coordinates {x, y, z) are adopted. Consider a quasi- 
neutral plasma with singly-charged ions of mass rrii. The ion/electron number density tiq is 
assumed to be uniform and constant. Suppose that Tj = rTf,., where Tj^e is the ion/electron 
temperature, and r is uniform and constant. Let there be no variation of quantities in the 
z-direction: i.e.., d/dz = 0. Finally, let all lengths be normalized to some convenient scale 
length a, all magnetic field-strengths to some convenient scale field-strength Ba, and all 
times to a/Va, where Va = Ba/ ^/JIonornl. 

We can write B = Vip x z + {Bq + bz) z, and P = Pq — B^bz + 0{1), where B is the 
magnetic field, and P the total plasma pressure. Here, we are assuming that Bq ^ Pq ^ 1, 
with ip and bz both 0(1).^^ Let, (3 = T Pq/Bq be (F times) the plasma beta calculated with 
the "guide- field" , Bq, where F = 5/3 is the plasma ratio of specific heats. Note that the 
above ordering scheme does not constrain f3 to be either much less than or much greater 
than unity. 

We adopt the reduced, 2-D, two-fluid, drift-MHD equations derived in the companion 



paper: ^^ 



^ = [<p-d,Z,^]+v{J-Jo) _ /fiMi±l) v2[\4 + {d,/c,)J], (1) 

Ot Cfj 

'dt ~ 



Z] + Cf, [Vz + idfs/cf,) J, ^lj]+c^DY + /ie df, W\U - d^Y), (2) 



^ = [0,f/] - ^ {V'[<P,Z] + [U,Z] + [r,0]} + [J,V^] +/i, V2(f/ + rf^rF) 

+/ieV2(f/-rf^F), (3) 

dV 

-^ = [0, Vz] + Of, [Z, ^] + /i, V'Vz + /ie V'[Vz + {dp/cp) J], (4) 

where D = r]{l - (3/2) [r/(l + r)]) + k/^, U = W, J = V^^, and Y = V^Z. Here, 
c^ = y^/3/(l + /3), dfs = Cfsdi/y/1 + r, Z = bz/c^ y/l + r, di = (mj/noeVo)^^V«! and 



[A, 5] = VA X VB-z. The guiding-center velocity is written: V = V0 x i + ^/T+TVz■ 



Furthermore, f] is the (uniform) plasma resistivity, /ije the (uniform) ion/electron viscosity, 
K the (uniform) plasma thermal conductivity, and Jo{x) (minus) the inductively maintained, 
equilibrium plasma current in the z-direction. The above equations contain both electron 
and ion diamagnetic effects, including the contribution of the anisotropic ion gyroviscous 
tensor, but neglect electron inertia. Our equations are "reduced" in the sense that they do 
not contain the compressible Alfven wave. However, they do contain the shear-Alfven wave, 
the magnetoacoustic wave, the whistler wave, and the kinetic-Alfven wave. 

B. Plasma equilibrium 

The plasma equilibrium satisfies d/dy = 0. Suppose that the plasma is bounded by 
rigid walls at a; = ±a;^, and that the region beyond the walls is a vacuum. The equilibrium 
magnetic fiux is written ^^^\x), where ^^'^\—x) = ijj^^\x), and (P%jj^^\x)/dx'^ = Jo{x). The 
scale magnetic field-strength, Ba, is chosen such that ^jj^^\x) -^ —x'^/2 as |a;| -^ 0. The 
equilibrium value of the field Z takes the form Z^'^\x) = —[V}y /dp (1 + r)] x, where V}y 
is the (uniform) total diamagnetic velocity in the y-direction. The equilibrium value of the 
guiding-center stream-function is written (p^'^^x) = —VJ^j^x, where VJ^^ is the (uniform) 
equilibrium E x B velocity in the y-direction. Finally, the equilibrium value of the field V^ 
is simply V"J°^ = 0. 

C. Asymptotic matching 

Consider a tearing perturbation which is periodic in the y-direction with periodicity 
length /. According to conventional analysis, the plasma is conveniently split into two 
regions. ^^ The "outer region" comprises most of the plasma, and is governed by the equations 
of linearized, ideal-MHD. On the other hand, the "inner region" is localized in the vicinity 
of the magnetic resonance x = (where M"^ = 0). Non- linear, dissipative, and drift-MHD 
effects all become important in the inner region. 



In the outer region, we can write ip{x, y, t) = ■i/j^^\x)+iIj^^\x, t) exp(i fcy), where k = 2ti /I 
and |^(i)| < |^(o)|. Linearized ideal-MHD yields [i)^^\ j(°)] + [^(°), J^^)] = 0, where J = V^^. 
It follows that 

The solution to the above equation must be asymptotically matched to the full, non-linear, 
dissipative, drift-MHD solution in the inner region. 

III. INTERACTION WITH A RESISTIVE WALL 
A. Introduction 

Suppose that the walls bounding the plasma at x = ±x^ are thin and resistive, with time- 
constant r^. We can define the perfect- wall tearing eigenf unction, ippw{x), as the continuous 
even (in x) solution to Eq. (5) which satisfies ippw{0) = 1, and ■?/'p^(±x^) = 0. Likewise, 
the no- wall tearing eigenf unction, ipnw{x), is the continuous even solution to Eq. (5) which 
satisfies ippwiO) = 1, and ■?/'p^(±oo) = 0. In general, both ippwi^x), and ipnw{x) have gradient 
discontinuities at x = 0. The quantity Apw = [dtjjpw / dx]Q'^ is the conventional tearing 
stability index^^ in the presence of a perfectly conducting wall {i.e., r^ -^ cxd), whereas 
^nw = [d'lpnw/dxJQll > Apw is the tearing stability index in the presence of no wall {i.e., 
T-u) -^ 0). Finally, the wall eigenf unction, ip^{x), is defined as the continuous even solution 
to Eq. (5) which satisfies ipw{^) = 0, ■ipw{±Xuj) = 1, and ■?/'^(±oo) = 0. This eigenfunction 
has additional gradient discontinuities at x = ±x^. The wall stability index, A^ < 0, is 
defined A^ = [di}^/dx]llt. 

According to standard analysis,'^ the effective tearing stability index. A' = [d\nip/dx]Qt, 
in the presence of a resistive wall is written 

T/2 /\ +T/2/\ 



where V is the phase- velocity of the tearing mode in the lab. frame, and Ku = {—'^w)/{k r^). 
Also, the net ^/-directed electromagnetic force acting on the inner region takes the form 

fy = -^{^n^-^P^.)^^^,^', (7) 

where ^{t) = \^^^\0,t)\ is the reconnected magnetic flux, which is assumed to have a very 
weak time dependence. 

B. Island geometry 

In the inner region, we can write 

ij{x,9,t) = + ^{t)cos9, (8) 

where 9 = ky. As is well-known, the above expression for ip describes a constant--?/' magnetic 
island of full- width (in the x-direction) W = Aw, where w = y^ . The region inside the 
magnetic separatrix corresponds to ■?/' > —^, whereas the region outside the separatrix 
corresponds to ^^ < —^. It is convenient to work in the island rest frame, in which d/dt ^ 0. 
It is helpful to define a flux-surface average operator: 

/ft I aw //(S'^'^) d9 , . 

(/(3,V^,^)) = /^^^ (9) 

for ip < —\I/, and 

J~eo 2 \x\ zvr 

for ip > —^. Here, s = sgn(x), and x{s,iIj,9o) = (with tt > ^q > 0)- The most important 
property of this operator is that {[A, ip]) =0, for any field A{s, ip, 9). 

C. Ordering scheme 

In the inner region, we adopt the following ordering of terms appearing in Eqs. (l)-(4): 



5 J = 5J^^\s,ij,9). Moreover, V = V(°), r = r^"), cp = cf , dp = df\ /i,,e = /ig, k = k^'^\ 
1] = r]^'^\ and d^/dt = d\P^'^ydt. Here, the superscript *^*^ indicated an ith order quantity. 
This ordering, which is completely self-consistent, implies weak {i.e., strongly sub-Alfvenic 
and sub-magnetoacoustic) diamagnetic flows, and very long {i.e., very much longer than the 
Alfven time) transport evolution time-scales. 

To lowest and next lowest orders, Eqs. (l)-(4) yield: 



;^ = [0(^) - d, Z('\^] + r^(^) 5J(^) - ^-'^^^^^ + "^ V^[V;(^) + {d,/c,) 6J(% 



^^ cos^= [(j)^''' -dpZ^^>,^\+r]'^'>6J^'> - ''' "^^" ' '' V'\Vy> + ida/ca)6J'^'>\, (111 
dt cp 

= c, [Vi^) + {d,/c,) <5J(2), ^] + c| D(^) r(i) + /i(2) rf, V^(^(i) - d, Y(% (12) 

+/ir V2(f/« + rf^rrW) + /ii^) V2(f/« - rf^F«), (13) 

= -M« [Vj'\^] + c, [Z(^), V^] + ^f ) VVi^) + /.f V^[V;(^) + {d,/c,) 5J(^)] (14) 

in the inner region, where (5J(2) = J + 1, F(i) = V^Z^^), f/^^) = VV^^\ M(^)(s,^) = 
dcp^^^dijj, and L*^^)(s,-?/') = dZ^^^/dil). Here, we have neglected the superscripts on zeroth 
order quantities, for the sake of clarity. In the following, we shall neglect all superscripts, 
except for those on 0*^^^ and Z^'^\ for ease of notation. 

D. Determination of flow profiles 

Flux surface averaging Eqs. (12) and (13), we obtain 

(V^f/) + '^P^^'^^-^'-) iy^Y) = 0, (15) 

(/ii + lie) 

and 

6^w^{V^Y)-{Y)=Q, (16) 

where 



6=^^.^^. (17) 

W VD V lii + fie 



In the following, we shall assume that 5^1. 

Now, we can write V^ ~ d"^ /dx^, provided that the island is "thin" (i.e., w <^ I). It 
follows that 

Mis, i,) = -MtlLZJ^ Lis, i,) + F(s, V^), (18) 



where 



dtp 



and 



^ (,...(,.) I) _(,.)i 



(f /, 4, dF\ 



(19) 



Note that Lis,ip) and F(s, ■?/') are ocfci functions of x. We immediately conclude that 
Lisjip) and Fisjip) are both zero inside the island separatrix (since it is impossible to 
have a non-zero odd flux-surface function in this region). The function L(s, ■?/') satisfies the 
additional boundary condition x L -^ V}y /d/s (1 + r) as \x\/w -^ oo. Here, we are assuming 
that w <^ Xyj. Moreover, the function F(s, ■?/') satisfies the additional boundary condition 
xF — > i\x\/xw) iV^'^^ —V) as \x\/w -^ 0, where V^'^^ is the unperturbed island phase- velocity 
in the lab. frame. 

It is helpful to define the following quantities: i(j = —ip/^, ((■ ■ ■)) = (■■■) w, and X = x/w. 
The solutions to Eqs. (19) and (20), subject to the above mentioned boundary conditions, 
are 

s\/i°^ 1 

^(^,^) = .... 7, ,^ 77V^> (21) 

and 



wdpil+r) ((X2))' 



respectively. Of course, both L(s,-?/') and F(s,-?/') are zero inside the island separatrix (i.e., 
ip < 1). In writing Eq. (21), we have neglected the thin boundary layer (width, 6w) which 



resolves the apparent discontinuity in L{s,iIj) across the island separatrix. This boundary 
layer, which need not be resolved in any of our calculations, is described in the companion 
paper. -"^^ Note that the function L{s,iIj) corresponds to a velocity profile which is localized in 
the vicinity of the island, whereas the function F{s,iIj) corresponds to a non-localized profile 
which extends over the whole plasma. 

E. Force balance 

The net electromagnetic force acting on the island region can be written^^ 

r — OD 

fy = -2k^ {SJs sine) dij, (23) 

where SJg is the component of SJ with the symmetry of sin 6. Now, it is easily demonstrated 
that 

{5J, sine) = ^{x[5Js,ij]), (24) 

so it follows from Eq. (13) that 

(„,.,,). _(if^_l((,»)0_2(.3>^-(.)f). (25) 



Hence, 



/ d^F dF ^ 

/, = 2(/i, + /ie) lim i{x') -2{x^) — -{x)F 

X/W-^QO \ dip dip J 



2s{j2i + /ie) lim 

x/w—>oo 



2 d /I d{xFy 
dx \x dx 



(26) 



Finally, Eq. (22) yields 






Equating Eqs. (7) and (27), we obtain the island force balance equation: 

^ '^' + '""""" -^^' ^^^.■-^^)^;^(.W. (28) 

Xm / V -r V,„ 



This equation describes the competition between the viscous restoring force (left-hand side) 
and the electromagnetic wall drag (right-hand side) acting on the island, and determines the 
island phase-velocity, V, as a function of the island width, W. Note that the above force 
balance equation is identical to that obtained from single-fluid MHD theory/ 

F. Determination of ion polarization correction 

It follows from Eqs. (11), (13), and (14) that 

"'--2 (^ - \lF j * I'-' "'' + "' ' ^" + " IF W- '''' 

where 5Jc is the component of 5 J with the symmetry of cos^. In writing the above expres- 
sion, we have neglected any boundary layers on the island separatrix, since these are either 
unimportant or need not be resolved in our calculations (see Ref. 16). Now, making use of 
Eqs. (18), (21) and (22), we can write 

M{s, ip) = '- C[i)) + J^(^), (30) 

and 

s(\/(o) -1/.^°^) . s(V^^)-V) 
M{s,ij) + dpTL{x,ij) = ^ ^£(V^) + ^^ —J'W. (31) 

W Xw 

Here, V^Jj^ = ( V?/ +''" K''!)'' )/(! + ''") i^ ^^^ unperturbed Ex B velocity, V^y the unperturbed 
ion velocity, and V^y^ the unperturbed electron velocity. [Note that V}y = V^y — Vj|J^.] 
Furthermore, V^^^ = {^iV^y + ^eV^^)/{^i + /ie) (see Ref. 16) is the unperturbed island 
phase-velocity, and V the actual phase- velocity. All of these velocities are measured in the 
lah. frame. Finally, both /^(V^) and J^{ip) are zero for ip < 1, whereas 

m = ^y (32) 



and 



10 



in the region ip > 1. 

Now 

4 r°° 
A\V) = - {{5Jc cos^))# (34) 

(see Ref. 14), where A'{V), which is specified in Eq. (6), is the effective tearing stabihty 
index in the presence of the resistive walL Hence, it follows from Eqs. (29), (30), (31), and 
(34) that 

(0) ^ n/(o) ym-^ 



7] dt ^ ^^ 2 ^^^^^3 



2(\/W-[\/i°i, + \//°V2)(\/W-\/) (\/W - Vf 



(35) 



where 



h = 2f^^^^^d^lJ = 0.823, (36) 

Equation (35) is the Rutherford island width evolution equation^^ for a propagating 
magnetic island interacting with a resistive wall. There are three separate ion polarization 
terms on the right-hand side of this equation. The first (second term on r.h.s.) is the drift- 
MHD polarization term for an isolated island (see Ref. 16), and is unaffected by wall braking. 
This term, which varies as W^^, is stabilizing provided that the unperturbed island phase- 
velocity lies between the unperturbed ion fluid velocity and the unperturbed E x B velocity, 
and is destabilizing otherwise. The third (fourth term on r.h.s.) is the single-fluid MHD 
polarization term due to the island velocity-shift induced by wall braking (see Ref. 9). This 
term is always destabilizing, and varies as W'^^ and the square of the wall-induced velocity- 
shift. The second (third term on r.h.s.) is a hybrid of the other two polarization terms. The 

11 



sign of this term depends on many factors. However, in the hmit of small electron viscosity 
(compared to the ion viscosity), when the unperturbed island phase- velocity lies close to the 
unperturbed velocity of the ion fluid, ^^ the hybrid term is stabilizing provided V}y V^^^ > 0, 
and destabilizing otherwise. In other words, the hybrid term is stabilizing if the unperturbed 
island propagates in the ion diamagnetic direction (in the lab. frame), and destabilizing if 
it propagates in the electron diamagnetic direction. The hybrid polarization term varies as 
W'"^, and is directly proportional to the wall-induced island velocity-shift. 

IV. INTERACTION WITH A STATIC EXTERNAL RESONANT MAGNETIC 

PERTURBATION 

A. Introduction 

Let the walls bounding the plasma at x = ±Xw now be non-conducting {i.e., r^ -^ 0). 
Suppose that an even (in x) static magnetic perturbation (with the same wave-length as the 
magnetic island in the plasma) is generated by currents flowing in fleld-coils located in the 
vacuum region beyond the walls. 

The no-wall tearing stability index, Anw, is defined in Sect. Ill A. The coil eigenfunction, 
ipc{x), is the continuous even solution to Eq. (5) which satisfies ipd^) = and ipd^Xw) = 1- 
In general, this eigenfunction has a gradient discontinuity at a; = 0. It is helpful to define 
A, = [d^ljJdx]'ot. 

According to standard analysis,'^ the effective tearing stability index. A' = [d\ntp/dx]Qll, 
in the presence of an external magnetic perturbation is 

A'{t) = An^ + ^c -^ COS^(t), (40) 

where ^{t) = \^^^\0,t)\ is the reconnected magnetic fiux, which is assumed to vary slowly 
in time, and ^c the flux at the walls solely due to currents flowing in the external coils. 
Furthermore, ip{t) is the phase of the island measured with respect to that of the external 
magnetic perturbation. Since the external perturbation is stationary, it follows that 

12 



|^*m (41) 

where V{t) is the instantaneous island phase- velocity. Also, the net y-directed electromag- 
netic force acting on the island takes the form 

fy{t) = --A,^^, sin ifit). (42) 

Note that, unlike the braking force due to a resistive wall, this force oscillates in sign as the 
island propagates. 



B. Determination of flow profiles 

We can reuse the analysis of Sect. HID, except that we must allow for time dependence 
of the function F to take into account the oscillating nature of the locking force exerted on 
the island by the external perturbation. Hence, we write 

Mis, ^, t) = _MtlLZf^ Li,^ ^) + Fis, ^, t), (43) 

[Hi + He) 



where 



sV, 



(0) 



wdp{l +r) ((X^)) 



and 



d_ 

dil> 



. , X <9 / . 4, dF\ . 2, dF 



0. (45) 



In order to proceed further, we adopt the separable form approach to solving Eq. (45) 
which was introduced and justified in Ref. 14. In other words, we try the following solution: 

F{s,ilj,t) = sFi{ilj) sin (f kV{t') dt') + s F2{ip) cos ( f kV{t')dt'Y (46) 

Of course, Fi[%Ij) and F2{iIj) are both zero within the island separatrix. Furthermore, 

klFi^Fo, (47) 

\x\ F2 -^ 0, (48) 
13 



as \x\/w -^ oo. Here, Fq is a constant. The above boundary conditions imply that the 
function F[s,ilj,t) corresponds to a velocity profile which is localized in the vicinity of the 
island. 

Matching to the outer region yields 

Fo sin n kV{t') dA = \/(°) - V{t). (49) 



Hence, differentiating with respect to t, we obtain 

1 dV 
kV 

and 



^ = - Fo cos (J k V{t') dA , (50) 



^(l\7^]=kV{V^'^-V). (51) 

dt\kV dt ) ^ ' ^ ' 

Substituting Eq. (46) into Eq. (45), and integrating once in i\} using the boundary con- 
ditions (47) and (48), we get 



2w'' dipy dip J 

2w^ dipy dip J w 



sgn(V) 2L JL f ((a-» ^"I - ((A-)> F,^-^. (53) 



Here, A = J2{iJLi + iie)/k\V\ is the localization scale-length of the velocity profile corre- 
sponding to the function F . 

Suppose that w <^ \ <^ Xyj. In other words, suppose that the localization scale- length of 
the velocity profile associated with F is much larger than the island width, but much smaller 
than the extent of the plasma. In this limit (which corresponds to the "weakly localized" 
regime of Ref. 14), Eqs. (52) and (53) can be solved to give 



w 



1-exp -^^ cos^ 



H^), (54) 



|X|F2 = sgn(\/)^exp('-^') sinj'^') J^{iP). (55) 

Here, J-'{ip) is specified in Eq. (33). It follows from Eqs. (46), (49), and (50) that 
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F{sJ,t) = -iV^^^ -V) 



1 — exp 



w\X\ 



A 



cos 



w\X\ 



\X\ 



s 1 dV 



w k \V\ dt 



exp 



w\X\ 
A 



sm 



'w\X\\ T{i)) 



A 



IXI 



(56) 



C. Island equation of motion 



Reusing the analysis of Sect. HIE, taking into account the time dependence of F, we 
obtain 



fy = 2s{iii + He) lim 

x/w~*oo 



X 



d (I d{xFY 



4, 



{x') ^ - (x) FJ #. (57) 



dx\x dx 

According to the boundary conditions (47) and (48), the first term on the right-hand side is 
identically zero. Transforming the second term on the right-hand side, using the fact that 
the integral is dominated by the region |X| ^ 1, we get 



d r^ diX F) 



(58) 



Finally, Eqs. (50), (51), and (56) yield 

'dV 



fy = ^ 



dt 



+ A;|\/|(\/-\/(°)) 



(59) 



Making use of Eq. (42), the island equation of motion takes the form: 



Here, (lVc/4)^ = ^c^c- The first term on the left-hand side represents the inertia of the 



region of the plasma (of width w2 (/Xj -|- iie)/k \V\) which is viscously coupled to the island, 
the second term represents the viscous restoring force, and the third term represents the 
locking force due to the external perturbation. Note that the above equation is identical 
to that obtain from single-fiuid MHD theory. ^^ The above analysis is valid provided w -C 

^2(/ii + IJ,e)/k\V\ < Xy,. 
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D. Determination of ion polarization correction 



Reusing the analysis of Sect. IIIF, we obtain 



where 



M{s,i^,t) = -^^^^^-^^^ C{^) - ;^f^ T(^), (62) 



and 



s 



a:!°i-o^M, ^/.(^) 



M(s,^,t) + rf^rL(x,^) = ^ ^£(^) - /''' ^ W- (63) 

Here, use has been made of Eqs. (56) and (59), as well as the fact that the polarization term 
integral is dominated by the region |X| ~ 0(1). Finally, Eqs. (34), (40), and (42) yield 



, A; (\/W - [\/i"i + <V2) fW^\' . 

U r — — Sm 09 

2 U^ + ^e) V 4 y ^ 



where Ji, /2, Is, and /4 are specified in Sect. IIIF. 

Equation (64) is the Rutherford island width evolution equation for a propagating island 
interacting with a static external resonant magnetic perturbation. There are three separate 
ion polarization terms on the right-hand side of this equation. The first (third term on r.h.s.) 
is the drift-MHD polarization term for an isolated island (see Ref. 16), and is unaffected by 
the external perturbation. The third (fifth term on r.h.s.) is the single-fluid MHD polariza- 
tion term due to the oscillation in island phase- velocity induced by the external perturbation 
(see Ref. 14). This term modulates as the island propagates, but is always destabilizing. The 
second (fourth term on r.h.s.) is a hybrid of the other two polarization terms. 
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E. Solution of island equations of motion 

Let us solve the island equations of motion, (41) and (60), in the limit in which the 
external magnetic perturbation is sufficiently weak that it does not significantly perturb the 
island phase-velocity. Let us also assume that rj is so small that the island width, W, does 
not vary appreciably with island phase. In this limit, we can write 

^(t) = k \/(°) t + as sin(A; 1^^°) t) + a^ cos{k \/(°) t), (65) 

where \as\, {ad ^ 1- Substitution of the above expression into Eqs. (41) and (60) yields 

f)7^)74A[K<»f, (66) 



and ftc — sgn{V^^^)as, where A = J2{fii + fXe)/k\V^'^d is the velocity localization scale- 
length. Averaging over island phase, using Eq. (65), we obtain 

COS^~-y, (67) 

sm^^sgn(V^(°))y, (68) 



sinV--- (69) 

Hence, the average of the Rutherford island width evolution equation, (64), over island phase 
takes the form 



hdw_ , , {v^'^-v^%){v^'^-vt ^ 



3 



The first two terms on the right-hand side of the above equation are the intrinsic tearing 
mode drive and the drift-MHD polarization term, respectively, and are unaffected by the ex- 
ternal perturbation. The next three terms (within the curly braces) are the phase-averaged 
external perturbation drive, hybrid polarization term, and single-fiuid MHD polarization 
term, respectively. It can be seen that the external perturbation drive is on average sta- 
bilizing, whereas the single-fiuid MHD polarization term is destabilizing^ The sign of the 
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hybrid term depends on many factors. However, in the hmit of small electron viscosity 
(compared to the ion viscosity), when the unperturbed island phase- velocity lies close to the 
unperturbed velocity of the ion fluid, ^^ the hybrid term is on average stabilizing provided 
V}y V^^^ > 0, and destabilizing otherwise. In other words, the hybrid term is stabilizing if 
the unperturbed island propagates in the ion diamagnetic direction (in the lab. frame), and 
destabilizing if it propagates in the electron diamagnetic direction. Finally, since our analysis 
is based on the fairly reasonable assumption that w/X <^ 1, it follows from Eq. (70) that 
the phase-averaged external perturbation drive dominates the phase-averaged hybrid and 
single-fluid MHD polarization terms. Hence, we conclude that, on average, an island prop- 
agating in the presence of an external magnetic perturbation experiences a net stabilizing 
effect. 

V. SUMMARY AND DISCUSSION 

We have investigated the dynamics of a propagating magnetic island interacting with a 
resistive wall or a static external resonant magnetic perturbation using two-fluid, drift-MHD 
theory in slab geometry. In both cases, we find that the island equation of motion takes 
exactly the same form as that predicted by single-fiuid MHD theory (see Sects. HIE and 
IV C). However, two-fiuid effects do give rise to additional ion polarization terms in the 
Rutherford island width evolution equation. 

In general, we find that there are three separate ion polarization terms in the Rutherford 
equation (see Sects. IIIF and IV D). The first is the drift-MHD polarization term for an iso- 
lated island, and is completely unaffected by interaction with a resistive wall or an external 
magnetic perturbation. Next, there is the polarization term due to interaction with a resis- 
tive wall or magnetic perturbation which is predicted by single-fluid MHD theory. This term 
is always destabilizing. Finally, there is a hybrid of the other two polarization terms. The 
sign of this term depends on many factors. However, in the limit of small electron viscosity 
(compared to the ion viscosity), when the unperturbed island phase- velocity lies close to 



the unperturbed velocity of the ion fluid/^ the hybrid term is stabihzing if the unperturbed 
island propagates in the ion diamagnetic direction (in the lab. frame), and destabilizing if it 
propagates in the electron diamagnetic direction. 

It is also demonstrated that a propagating magnetic island interacting with a static 
external resonant magnetic perturbation generally experiences a net stabilizing effect (see 
Sect. IV E). This follows because in the Rutherford island width evolution equation the 
phase-averaged drive term due to the external perturbation (which is stabilizing) is gen- 
erally much larger than either the phase-averaged hybrid polarization term (which can be 
destabilizing) or the phase-averaged single-fluid MHD polarization term (which is destabi- 
lizing) . 
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